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ABSTRACT 

I solve analytically the viscous evolution of a self-irradiated accretion disc, as seen 
during outbursts of soft X-ray transients. The solutions predict steep power-law X- 
ray decays Lx ~ (1 + V^visc)"'*, changing to Lx ~ (1 — i/^visc)'^ ^^^'^ times, where 
i vise, i vise ^-re viscous timescales. These forms closely resemble the approximate expo- 
nential and linear decays inferred by King and Ritter (1997) in these two regimes. The 
decays are much steeper than for unirradiated discs because the viscosity is a function 
of the central accretion rate rather than of local conditions in the disc. 
Key Words: accretion, accretion discs - instabilities - stars: X-rays: stars. 



1 INTRODUCTION 

Soft X-ray transients (SXTs) are close binary systems with 
, black-hole or neutron-star primaries; they undergo large 
outbursts, reaching X-ray luminosities of order the Edding- 
ton limit. The outbursts have some resemblance to those of 
dwarf novae, in which the primary is a white dwarf, with 

■ the important difference that the timescales are far longer. 

, SXT outbursts last of order a year, and recur at intervals of 

■ several years or greater, whereas these timescales are typi- 
' cally days and weeks for dwarf novae. It is by now generally 
, accepted that dwarf nova outbursts result from a thermal- 
' viscous disc instability (see e.g. Cannizzo 1993 for a review). 

While this picture naturally explains observed dwarf nova 
, timescales it requires essentially ad hoc modification to ex- 
• plain SXT light curves. However a crucial observed differ- 
] ence between dwarf nova disc and those in SXTs is that 

■ the latter are heavily irradiated by the central X-ray source 
' (van Paradijs & McClintock, 1994; Shahbaz & Kuulkers, 
' 1997). A recent paper (King & Ritter, 1997, hereafter KR) 

has pointed out that this fact inevitably lengthens the entire 
SXT outburst cycle, as irradiation will prevent the encroach- 
ment of the cooling front which terminates a dwarf nova 
outburst. Instead, the disc can only return to the cool state 
as its central accretion rate Mc is reduced, which in turn 
requires most of the disc mass to be accreted. KR modelled 
this effect by simple dimensional relations and showed that 
Mc (and thus the X-ray emission) will decay roughly expo- 
nentially on the hot-state viscous time of the irradiated disc. 
This prediction, and the related timescales and luminosities, 
are in good agreement with observed SXT light curves (see 
Tanaka & Shibazaki, 1996 and Chen, Shrader & Livio, 1997 
for recent reviews). 

KR pointed out that long decays were a feature of the 
now-discarded mass transfer burst model of dwarf novae, 
which shared the property that the outburst could only 
terminate through the accretion of most of the disc mass. 
However Cannizzo (1997) has argued that since the mass 
transfer burst model gives rather shallow power-law decays 



(typically Mc ~ with e << 1; Pringle, 1974; Can- 

nizzo, Lee & Goodman, 1990) rather than exponentials, the 
irradiated-disc model of KR would give a poor representa- 
tion of observed SXT light curves. I consider this question 
further in this paper. 

In the unirradiated disc models considered by Pringle 
(1974) and Cannizzo, Lee & Goodman, (1990) the viscos- 
ity controlling the time evolution of the disc is a function 
of purely local variables, either through assumed power-law 
dependences (Pringle) or the alpha-prescription 



f = acsH 

1/2 

with the sound speed cs oc and scaleheight 



H 



cs 



\GM 



1/2 



(1) 



(2) 



fixed by the disc midplane temperature Tm resulting from 
local viscous dissipation (here R is the radial disc coordi- 
nate and M the central mass) . The irradiated-disc model of 
KR is fundamentally different; the temperature and viscos- 
ity are no longer fixed locally, but instead by Mc, and have 
explicit time dependences. I shall find analytic solutions of 
the disc diffusion equation for this case. These give much 
steeper power-law behaviour for Mc, which is in practice 
indistinguishable from the exponential decay predicted by 
KR. 



2 IRRADIATED DISCS 

The surface temperature of a disc irradiated by a central 
source is given by 



?7Mcc2(l-/?) 




'cUnH 


47rcri?2 


yRj 


dlni? 



(3) 



where r; is the efficiency of rest-mass energy conversion into 
X-ray heating, Mc the central accretion rate, H the disc 
scaleheight at disc radius R, (3 the albedo of the disc faces. 
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and the factor in square brackets lies between 1/8 and 2/7. 
The index n = 1,2 for irradiation by a central point source or 
the inner disc respectively (van Paradijs, 1996; Fukue, 1992; 
King, Kolb & Szuszkiewicz 1997). Neutron-star accretion 
always has n = 1, while n = 2 is possible in the black 
hole case because of the lack of a hard stellar surface; this 
effect probably allows black-hole binaries to be transient at 
mass transfer rates which would otherwise give a persistent 
system (King, Kolb & Szuszkiewicz 1997). However during 
outbursts, observations of a strong power-law continuum 
strongly suggest that the disc develops an extended X-ray 
corona, thereby reverting to n = 1. The distinction between 
these two possibilities will in any case be unimportant for 
the purposes of this paper (I shall assume n = 1 below): 
the main effect comes from the radial dependence of Tin-, 
which is the same regardless of n. The ratio H/R is roughly 
constant in a disc, so Tin falls off as R~^^^. Thus for a disc 
with a large enough ratio of outer radius Rq to inner radius 
Rt, Tirr dominates the disc's own effective temperature T^s, 
given by 



4 _ 3GMM 



(4) 



(e.g. Frank et al., 1992) which goes as R-''^^'^ (here M is the 
accretion rate at radius R). Thus 



7^ = 3:r:^^(^-^):r 



A\nH 
A\nR 



(5) 



where Rs = 2GM/(? is the Schwarzschild radius of the cen- 
tral object. In X-ray binaries one typically finds H/R ~ 
0.2, /3 = 0.9 (de Jong, van Paradijs & Augusteijn; 1996), 
while the last factor on the rhs lies between 2/7 and 1/8. 
Further, the inner disc radius R^ is close to Rs in all X-ray 
binaries, so the efficiency is JiJO.l. With a typical value 
T] = 0.2 we find 



r.4 

irr 



6.7x10-"^ 
Rs 



(6) 



for a quasi-steady disc {M ~ Mc). Since we typically have 
Rs^2 X l(f cm and a disc radius _Ro > 10^^ cm, we find 

r4/r*fl>34. 

This result shows that irradiation dominates the output 
of the disc over most of its surface in any situation where 
the central luminosity is significant. In fact we can go further 
and show that irradiation completely alters the disc struc- 
ture over most of its area. For Tirr generally exceeds the 
midplane temperature Tm of the outer disc, given by radia- 
tive equilibrium as ~ tT^s ~ lO^^lf . Here r ~ 10 is the 
vertical optical depth of the unirradiated (cool) outer disc, 
with typical radu Ro > lO" cm (e.g. Frank et al., 1992). The 
outer disc therefore develops an almost isothermal vertical 
structure, and we can replace Tm in eg by 

Tirr « M^^^R-^^^ 



SO that (\l 



1/4 



(7) 



(8) 



where z/* , M* are the viscosity and accretion rate at the inner 
disc radius _R, at time t = 0, measured from the beginning 



of the decay (the X-ray peak). Then the diffusion equation 
determining the disc surface density E(_R,t) 



as 3 a / 1/2 d . r,i/2vi\ 



(e.g. Frank et al., 1992) can be rewritten as 

1/4 

dR 



(9) 



(10) 



Now defining p = 2i?^/^ and S = i?^/^E we get finally 



dS_ 

dt 



R, 



Mc_ 



Thus we obtain a diffusion equation for S{R,t), which is 
however nonlinear because of the factor Mc^^ on the rhs: 
this itself depends on S and thus t through the relation 

Mc ^ 3TTv{R,,t)j:{R,,t) (12) 

(e.g. Frank et al., 1992). Using (^) and the definition of S 
we have 



Mc 



371!/, 



Mc 



1/4 



R. 



-3/2 



S{R*,t),] 



(13) 



and thus 

Mc" 
M, 



3Tvu,S{R,,t) 



-I 1/3 



(14) 



It is clear that ( |ll| ) has a steady solution of the form S 

ap + b, or 

b 



a 

R ^ 



(15) 



Since from (g) u (x R, this is simply the usual steady solution 
i/E = A + BR'^''^ (e.g. Frank et al., 1992), in which the first 
term represents a steady inflow with accretion rate given by 
A, and the second term gives zero net mass inflow. 

I now seek separated time-dependent solutions of (|ll]) 
of the form 

s{R,t) ^ R:l'^e{t)P{p) (16) 

with 6 it) dimensionless and P with the dimensions of E. Eq. 
(0) gives 





'3ttu,P{p*)' 







1/3 



nl/3 



with p, = 2i?y^ Choosing 61(0) = l,P(p,) 



(17) 

E(i?.,0), the 



deflnition of A/, implies that the square bracket = 1, so that 

1/4 

^e^'^'it). (18) 



Then inserting (0) into (|ll|) and using disl ) gives 



R* P 



1 



(19) 



3614/3 

where T is a separation constant with dimensions of time. 
Since 6(0) = 1 we find 

-3 



1 + 



(20) 
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and also 



(21) 



Replacing p by 2R^^^ and assuming that the disc is irradi- 
ated out to a fixed outer edge Rq, i.e. 



E(i?o,t) = 0, 

we get the eigenvalues 



T 



and 



i^,(n- 1/2)2 7r2' 



n= 1,2,3, . 



(22) 



(23) 



(24) 



for the cosine and sine terms respectively. We can simplify 
these by noting from (H) that 



vi = v{Ro,0) = T^*^, 



so that 
T„ 

T' 



(2n- 1)2' 

^visc 



(25) 

(26) 
(27) 



(2n)2 

for the two cases, with ivisc = WRq/tv^pi a measure of the 
viscous time at Rq. Then in terms of the variable 

1/2 



TV ( R 



2 V Ro 



(28) 



we finally get the general solution of the diffusion equation 
for a fully irradiated disc obeying (^^ as: 



E(i?,t) = 

1 I A,i cos(2n — l)x 



B„ sin 2nx , 

+ r. . Z so,,, (29) 



a:3 ^ 1 [1 + (2n - l)2t/tvisc]^ [1 + {2n)H/Ui,c]'' 



Here the coefficients A„ , _B„ are to be found by ex- 
panding the initial surface density distribution E(_R, 0) as 
a Fourier series in cos(2n — l)x and sin2na;, viz 



4 r^^ 

- x^T,{R, 0) cos(2n ~ l)2;da; 

^ Jo 



Bn — — 



tt/2 



X E(_R, 0) sin 2na;da:. 



(30) 



(31) 



(If we had not used the boundary condition (^2|) we should 
have found a Fourier transform rather than a series.) For 
an initial density distribution with no discontinuities we 
therefore have A„, B„ ~ 1/n^. Note that the steady so- 
lution (^sj) does not obey the boundary condition (p^, 
as mass has to be fed in from R > Rt, to maintain the 
steady state. Accordingly if we start with initial conditions 
E(7?, 0) — ax~^ + bx~^ in ( |29| ) we do not obtain the steady 
solution, since the expansion assumes that ( p^ holds. It 
therefore demands a discontinuity in E(_R, 0) at i? = i?o (i.e. 
An,B„ ~ 1/n) which subsequently diffuses away. 

Given the solution (^9|) , the central accretion rate Mc (t) 
can be found from ([l7|). (I do not consider here the possible 
glitch in the light curve caused by irradiation of the outer 



disc (cf KR).) The disc mass within a radius given by xr at 
any time can be expressed as 



M(R,t) 



64i?g 



X E(a::,t)da;. 



(32) 



Using ( |2S| ) with x = 7r/2, the total mass in the nth eigen- 
mode is 



64i?g(-l)"+M„ 



+ 



64fig((~l)"+^ + l)g„ 
(2n)7r3[l-f (2n)2f/tviBcl = 



(33) 



The structure of the expansion (^9|) is simple to under- 
stand. The viscous time of the nth disc eigenmode varies 
as tn ~ tvisc/n^, reflecting the fact that steeper spatial 
gradients cause more rapid diffusion. While the full spec- 
trum of eigenmodes is in general present in the surface 
density distribution at the beginning of the outburst, the 
higher eigenmodes will decay more rapidly (as {l + t/t„)~'^) 
as the outburst proceeds, causing E to approach the flrst 
term of the expansion (cx cos a;). The disc mass is dom- 
inated by this mode in any case, as from ( p3[ ) we have 



E(7?,t) ~ Ai 



1 + 



Thus 

-3 



3/2 



1/2 



For this term, ( |17[ ) gives simply 

-4 



1 + 



t 



(34) 



(35) 



Accordingly the X-ray luminosity must decay at least as 
steeply as this. From (^5|), the typical 'exponential' decline 
phase of SXT light curves, during which the flux drops by a 
factor ~ 10, corresponds to t running from to 1.26iio/^i. 
The paper by KR derived the simple law Mc (x e"*''^ for the 
typical 'exponential' decline phase of SXT light curves, with 
r = Rg/Siyi. The best-fitting function of the form (BSl) to 
this exponential, over a dynamic range of 10'', has normal- 
ization 1.05 and tvisc ~ 3.11t = 1.035i?o/i^i, as opposed to 
ivisc ~ IGRii/'K^vi = l.G2R^/vi. Thus the simple treatment 
of KR agrees surprisingly well with the present detailed cal- 
culation. As can be seen from Fig. 1, the power-law decay 
is almost indistinguishable from the exponential over this 
range. It is much steeper than the Mc ~ f"'-^^''' decay of an 
unirradiated disc (see Section 4). 

I note finall y in this Section that the fundamental disc 
outburst mode (pi) has the mass concentrated quite close 
to the centre: from (B2h we find 



R_ 



1/2 



(36) 



Thus one-half of the mass of this mode (i.e. effectively one- 
half of the total disc mass) lies within R = 0.11-Ro- This 
contrasts with outbursts of unirradiated discs, where the 
usual E ~ i?"^/" distribution leads to M{R,t) ~ R^^^, so 
that much of the mass is still in the outer parts of the disc. 
The reason for the difference is the larger viscosity in the 
outer parts of an irradiated disc (see Section 4 below) , which 
drives the mass inwards more rapidly. 

This difference in mass distribution holds only in out- 
burst. In quiescence irradiation is of course completely neg- 



© 0000 RAS, MNRAS 000, 000-000 



4 A. R. King 



ligible. The disc will evolve in the usual way to a state in 
which most of the mass is near its outer edge (i.e. E ~ J?) 
and begin the next outburst in this state (cf Cannizzo, 1993). 



3 THE LATE DECLINE 

The solutions derived above assume that the whole disc is 
dominated by irradiation. However as explained by KR, at 
late times this cannot be true, and the region where irradia- 
tion is important will retreat inwards as, R (x Mc . Indeed 
sufficiently large discs will already be in this regime at the 
start of the outburst, as even the Eddington limit luminos- 
ity will be unable to ionize hydrogen out to the edge of such 
a disc. KR showed that the decline of Mc changes from ap- 
proximately exponential to linear in this regime. 

The variables (R, t) or {p, t) used in the diffusion equa- 
tion up to now are not useful in this regime, as one wants 
to set boundary conditions where the disc structure changes 
from irradiated to effectively unirradiated, i.e. at a time- 
dependent radius R cx Mc . Accordingly we change the 
variables from (p,t) to {u,t), with 



P 



Ml'\t) 

After some manipulation, the diffusion equation 
then be cast in the form 



dS_ _ uMc dS 



R,Ah 



dS 



+ 



(37) 



(38) 



One can again look for separated solutions, this time of the 
form 



which implies 

Mc 



uU' 



ml^^ U 



+ ■ 



3u, 



1 tr ir_ 

uU ^ u 



(39) 



(40) 



Evidently the only hope of finding separated solutions is to 
require 



Mc 



AM, 



3/4 



= A 



with A a constant. Then 
Mc = {At + Bf 

with B also constant. Substituting in 



gives 



{At + b)^=Au— + 

u r,mI'* 



iir_ 

uU ^ 



(41) 



(42) 



(43) 



which requires each side to equal a separation constant, say 
(3. Thus we find 



'^M^t + B) 



(44) 



Now (|l^ and (^IWequire (j) cx Mc'^ oc {At + Bf, which is 
compatible with~(fl4|) only if 13 — 3A. The u - equation can 
be rearranged as 

U" + -U' + CuU' - 3C(/ = 0, (45) 
u 

where C = AR^mI'''^ /Su,. This has a regular singular point 



at u = and an irregular singular point at u = oo, suggest- 
ing a relation to the confluent hypergeometric equation. A 
change of independent variable shows that 



2 



(46) 



where M{a, 6 : u) is the confluent hypergeometric (Kummer) 
function. The full solution is 

E(«,t)= j J{C){Ct + DfM(^-^,l-~^^AC, (47) 

recalling that C is proportional to A. We can use the freedom 
in /(C) to satisfy the boundary condition specifying Y,{uQ,t) 
on the irradiation boundary u = wo where Tirr = Tk. 

Ahhough (||) appears complex, its meaning is again 
straightforward, and similar to that of the expansion (E9|). 
The time development Mc oc {At + Bf shows that yl^ is 
a viscous timescale. Since Mc is decaying, we have approxi- 
mately 



Mc 



tL 



(48) 



where f^isc is the longest viscous time of the partially- 
irradiated disc. This expression predicts that the outburst 
of the irradiated disc shuts off completely at t = t^isc- In 
reality the unirradiated disc outside the irradiation bound- 
ary u = uo may produce further outburst activity. Note also 
that the innermost part of the disc may not be irradiation- 
dominated, as for small R the effective temperature T^s is 
larger than Tin- for similar local accretion rates. 

Equation (^) thus gives an approximate description of 
the regime in which the central emission is unable to keep the 
disc ionized all the way to its outer edge. KR found a linear 
decay of the central X-rays in this case, which is approxi- 
mately what is observed (Giles et al., 1996) in the decline 
of GRO J1744-28, the transient with the longest known or- 
bital period. As can be seen (Fig. 2), eq. ( p^ gives a rather 
similar light curve over the typical observed dynamic range. 
Note in particular that the linear decay of GRO J 1744-28 
extends only over a dynamic range of about 4 in flux. After 
this the observed fiux is persistently higher than the linear 
trend, just as in Fig. 2. 



4 DISCUSSION 

I have shown that analytic solutions for the viscous evolution 
of an irradiated disc predict a steep power-law decay in the 
central accretion rate Mc- In practice this is quite close to the 
exponential decay predicted by the dimensional arguments 
of KR; this is a good representation of typical observed X- 
ray light curves. It is important to ask why these solutions 
are so different from the very shallow Mc ~ t"''^^"' decay of 
an unirradiated disc. 

As pointed out by Cannizzo (1997), the latter is easy 
to understand. The equation for vertical radiative energy 
transport within a disc 



3kE 



9 ^GM 



(49) 



(where Tm is the midplane temperature) can be combined 
with the alpha-prescription (M to give 
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,1/3-02/3 



R, (50) 

where k is the opacity. The factor k^^^ generally varies only 
slowly with E and Tm, so the difltusion equation (^ predicts 



as 

dt 



oc E 



5/3 



(51) 



and thus E ~ t~^^^ for a separated solution for an unirradi- 
ated disc. 

The fundamental difference for an irradiated disc is 
that, from the viscosity u is independent of the local 
value of E. Thus (^l|) is replaced by 



^ocMy^E. 
dt 



(52) 



Were it not for the weak M^"^ dependence we would get 
an exactly exponential decay for E(t) and thus for Lx'- this 
is in any case clearly a good approximation until Mc has 
decreased by large factors. Note that this dependence of u 
on disc properties is global, i.e. on the value of M or E at the 
inner edge of the disc, rather than a local dependence on E 
at each radius. The effect is to convert the exponentials into 
steep power laws since the relation Mc oc E(_R*)^''^ implies 
Mc oc Mc^^ and thus Mc ~ for separated solutions. 

The different viscosity behaviour in an irradiated disc is 
also responsible for the greater central concentration of the 
disc mass noted at the end of Section 2. The unirradiated 
viscosity ( |50| ) drops as E decreases, tending to slow the de- 
pletion of the lower surface density of the outer disc regions. 
By contrast the irradiated viscosity ^ is independent of 
the local value of E and increases with R, thus driving mass 
inwards from large radii. The depleted density in the outer 
regions implies a smaller local accretion rate there: from (^ 
and the local equivalent of (^2|) we find 



M{R,t)^Mcit){ ^ 



1/2 



1/2 



cos f {R/Ro 
cosf (7?,/7?o)i/2' 



(53) 



so that the local effective temperature goes as 
stead of the usual ^ R~^^*. 



REFERENCES 

Cannizzo, J.K., 1993, in Accretion Disks in Compact Stellar Sys- 
tems, ed. J.C. Wheeler, World Sci. Pub. Co., Singapore, p. 6 

Cannizzo, J.K., 1997, to appear in Proceedings of the 13th North 
American Workshop on Cataclysmic Variables and Related 
Objects 

Cannizzo, J.K., Lee, H.M & Goodman, J., 1990, ApJ 351, 38 
Chen, W., Shrader, C.R. & Livio, M., 1997, ApJS, in press 
de Jong, J. A., van Paradijs, J, & Augusteijn, T., 1996, A&A 314, 
484 

Frank J., King A.R., Raine D.J., 1992, Accretion Power in Astro- 
physics, 2nd edn. Cambridge Univ. Press, Cambridge 
Fukue, J., 1992, PASJ, 44, 663 

Giles, A.B., Swank, J.H., Jahoda, K. Zhang, W., Strohmayer, T., 

Stark, M.J. & Morgan, E.H., 1996, ApJ 469, L25 
King, A.R., Kolb, U., & Burderi, L., 1996, ApJ 464, L127 
King, A.R., Kolb, U., & Szuszkiewicz, E., 1997, ApJ, in press 
King, A.R., Ritter, H., MNRAS, in press (KR) 
Pringle, J.E., 1974, Ph.D. thesis 
Shahbaz, T., Kuulkers, E., 1997, MNRAS, in press 
Tanaka, Y., Shibazaki, N., 1996, ARAA 34 607 
van Paradijs, J., 1996, ApJ, 464, L139 
van Paradijs, J. & McCUntock, J.E. 1994, A&A, 290, 133 

FIGURE CAPTIONS 

Figure 1. Comparison of the exponential light curve oc 
e~^^^ (solid curve) derived by KR, and the best-fitting 
curve of the form (^5|) of this paper (dashed curve): this 
has oc 1.05(1 + t/(3.11)r)-''. 

Figure 2. Comparison of the linear light curve oc l — t/r' 
derived by KR (solid curve), and the best-fitting curve 
of the form ( ^ ) of this paper (dashed curve): this has 
oc 1.073(1 -t/2.647r')^. 



5 CONCLUSIONS 

The presence of a central irradiating point source completely 
changes the viscous evolution of an accretion disc. Because 
the viscosity is controlled by conditions in the centre, the 
disc mass in outburst is more centrally condensed than for 
unirradiated discs. Most importantly, the decay of the cen- 
tral accretion rate is much steeper, being approximately ex- 
ponential. I conclude that irradiated discs give a good rep- 
resentation of soft X-ray transient outbursts. 
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